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1. Introduction 

Quantization of simple integrable systems is important as a test of quantization 
methods in more general cases. Any polynomial field operator must be well 
defined, and so must be regularized, respecting the demand of integrability. 
In ordinary integrable theory integrals of motion are polynomial field operators, 
and so have such infinities. Commutators of those integrals of motion also 
have their own infinities. Our problem is to determine those operators and 
a way to work with polynomial field operators which is free from all these 
infinities. We suppose that know the answer for quantum commutators such 
integrals of motion. If coupling constants don't take special values integrals 
of motion must commute. So our problem to construct such operators. The 
sine-Gordon and Zhiber-Shabat models are very convenient for investigating 
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this problem. 

In [1], two different methods of quantization are distinguished. The first 
is to construct the monodromy matrix T(A) and impose 

[T(\),T(jm)] = 0. 

The second is to construct the quantum integrals of motion Jj and demand 

[/,,/,] = <>. 

Classically, the methods are equivalent. At the quantum level, however, 
that is not the case. On the one hand, the generating function of quantum 
integrals of motion T(A) is well defined. On the other hand lnT(A) is not 
well defined, because we have a problem connected with divergences [1]. 
Let us note that for factorization of the S matrix, we must have quantum 
integrals of motion (which commute with the Hamiltonian), i.e. the second 
method of quantization is prefereed. This method have successfully tasted in 
quantum integrable sine-Gordon, Zhiber-Shabat system, quantum nonlinear 
Scrodinger equation, Korteweg-de Vries and modified Korteweg de Vries 
equations [1]. 

We have used this second method on the sine-Gordon and Zhiber-Shabat 
models and have obtained quantum integrals of motion for arbitrary coupling 
constants. If the coupling constants take special values, an interesting situation 
arises. Local conservation currents have been found and those currents have 
unusual commutation relations, i.e. new algebras. The generators of the 
algebras have two indices. This contrasts with ordinary conformal symmetry, 
where we have the Virasoro algebra with one- index generators [2]. Using this 
new symmetry, we can calculate correlation functions. 

The power law behavior of the resulting correlation functions is connected 
with the violation of the standard equivalence [3] between the sine-Gordon 
and massive Thirring fermionic models. Let us recall basic features of this 
standard equivalence, and explain probably hidden difficulties of such approach. 

We can consider work about equivalence of sine-Gordon/Thirring models 
[3] (work of Mandelstam and in perturbative approach works of Klaiber+Coleman) 
and explain our position. In our approach we use only Hamiltonian quantization. 
In the works of Mandelstam and Klaiber, another method has been followed. 
They solve the normal-ordered Lagrange equations and then check the equal- 
time anticommutators for fermionic fields. Reading those works carefully 
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reveals an equal-time anticommutator of the non-canonical form 

$(x, t)J + (y, t)] + = (x- yYS(x -y)± 5{x - y). 

We use the Hamiltonian approach and postulate the canonical equal-time 
anticommutators (for initial operators) at the beginning of our calculations. 

The Hamiltonian approach guarantees the conservation this equal-time 
anticommutator. It is easy to see this from the solution of the Hamiltonian 
equations 

ip(x, t) = exp {^—itH^j i/j(x, 0) exp (itH^J , 

where H is the quantum Hamiltonian of the system. So we don't have 
similar problems. Calculation of integrals of motion demands correct equal- 
time anticommutators because uses in fact basis property of 5 - function 
/ f(x)S(x — y)dx = f(y), and so we can not use above anticommutators as 
in the works of Mandelstam or Klaiber. 

In our opinion, the equivalence of these theories must be investigated in 
the Hamiltonian approach and a connection between integrals of motion for 
sine-Gordon and massive Thirring models must be found. 



2. Two Hamiltonian structures for the sine-Gordon equation 

It is well known that the classical model obeying the sine-Gordon equation 

d^(x,t) - d 2 x Mx,t) + (m 2 /[3) sm(fa(x,t)) = (2.1) 

is an integrable Hamiltonian system [4]. Let us introduce equal-time canonical 
Poisson brackets for (p(x,t),ir(x,t), where ir(x,t) = d t ip(x,t): 

{if(x,t),7r(y,t)} = 5(x-y). (2.2) 

Equation (2.1) can be written in Hamiltonian form with the Hamiltonian 

/+oo /II \ 

^ dx (-7T 2 + -{d xV f + (mV^Xl - cos(/V))J • 

Integrals of motion are determined in terms of (p(x,t),ir(x,t) by demanding 
that they commute with H', using the brackets (2.2). 
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The equation of motion also has a second Hamiltonian structure. Let 
us introduce the coordinates £ = x + t, rj — x — t, and define xiCiV) = 
</?((£ + 77) /2, (£ — rf)/2). In the new coordinates (2.1) becomes 

dMi, V) ~ (m 2 /P) sm(Px(Z, v)) = 0. (2.3) 
The canonical Poisson brackets for 77) are 

{x(e^),x(e'^)}=5(£-£') 
Equation (2.3) can therefore be written in Hamiltonian form with 

/+00 
^(cos(/3 X (£,77))-l). (2.5) 
-co 

The Hamiltonian and integrals of motion can be determined at some initial 
time. The condition on the initial function is a holomorphic condition d v x{!;, i]) = 
(77 is the time) if we have t = it. We should determine a quantum field 
of the theory as x(£,v) — exp(—ir}H) x(0 exp(ir]H). Integrals of motion 
of the Hamiltonian (2.5) and brackets (2.4) can be written like differential 
polynomials of <9gx(0- We will consider the following boundary condition 
d^xiO (£ ~~ *■ i 00 )) an d it has some singularity in a restricted region of the 
real axis £. Let us consider £ as a complex coordinate and expand d^x(0 m 
some series. Indeed, any holomorphic function in the region r < |£ — a\ < R 
can be represented by a series in positive and negative powers of (£ — a), the 
value of a-places where we have a singularity [5]. So we have 

and 

x (0 = 9+plne-E— ■ (2-6) 

We have set a = without loss of generality. Now let us consider quantization 
of our field. The commutation relations for the operators are 

\p, q] = 1, [a n , a m ] = n5 n+mfi - 

The normal ordering of non-commuting operators obeys 

pq =: pq : +1, qp =: pq :, 



(2.4) 
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Qj n (l— n . (X n (l— n • ~\~Ti^ CL— n Cl n . Cl— n CL n 71 ^> 0. 

The relation for field is the following 

Quantum integrals of motion are differential polynomials of x with some 
coefficients. These coefficients should be defined by 

[H,l k ]=0, (2.7) 

where the quantum Hamiltonian H of the theory is 

/+oo r+oo 
■ exp(«x(0) : + / : ex P : ^- 

-oo J —oo 

When we integrate over £, we obtain a singularity. Indeed, when we expand 
the product of ordered operators, we have, for example 

: exp(a*(0) = = ■= £ TT^W = exp(a*(0) ■ ■ 

We should determine the function at the singularity. Let us consider the 
commutator: 

[HMC)}=ni k (t)-ik(t) H, 

where ik(£) is the density of the integrals of motion. In the first piece we 
take the contour of integration which is above and in the second piece we 
take the contour of integration which is below Then we add integration 
over (Tr^oo which is equal to zero. The closed contour can be deformed and 
we obtain 

3. Quantum symmetry of integrable theory 



Using the definition of the commutator, we can calculate the integrals of 
motion. If the coupling constants take generic values, we have 

i A = : (dx) ■ + 2 : ( 5 X) : > 
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it = ■ Wi)' : ^'"^f ' 4 ' : m 3 ?X : + 

L55a' 
Ik* 1 " 



(3a 8 - 40a 6 + 155a 4 - 160a 2 + 48) _ 3 2 _ 



and they satisfy 

[J 4 ,J 6 ] =0. 

But if the coupling constant takes special values, some interesting results 
appear. We have obtained 

[n, J k m = o, 

if the coupling constant a = 1 in the sine-Gordon theory. Here we found 
J\0 = ^y(: (^ 1+fc expx(0)^exp-x(0) : 



+ (fc + 2 )(fc + i) : {d * +k 6XP m) 6XP : ' k E N - 

Similar cases are listed in the Appendix. 

We can calculate the commutation relations between the J% 

k tP~\~1' k—V 

[Jt J P J =(* + !)! £ y !(fc +T-r)! + ^ + 1 " /)_ 

-(p + i)iE^^n o (^^i-.) 

(p+l)!(fc + l)! 
+<W,o 2{k + p + 3)[ x 

(2+fe+p 2+fc+p \ 

n (n+k+i-j)-(-ir +i n (m+ P +i-j) , 
3=0 3=0 ) 

where = ^ 4b ^ fc (OC fc+1+n ^- This is a new type of infinite-dimensional 
algebra, with central extension. Using this quantum symmetry we can calculate 
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the correlation functions in a system of this sort. 



4. Evolution of operators in the Heisenberg picture 

In the Heisenberg picture, the operator 77) can be obtained from the 
operator A(£): 

A(^rj) = exp(-triH) exp^W). 
The correlation function of the operators has the form 

(0\A 1 ^ 1 ,r ]l )---A n ^ n ,r ]n )\0) = 

(0| exp(-i77iW) Ai(^i) expii^H) ■ ■ ■ exp(-ir} n H) A n (£ n ) exp(irj n H)\0) , 

where |0) is the vacuum of the theory. If 77) is a complicated function of 
operator x, we should fix a certain ordering. We choose, for example, 

exp(ax(£,??)) = exp(-irjii) : exp(«x(0) : exp^W), 

x{0 is the field of theory and it has the decomposition (2.6). The normal 
ordering : : was discussed in section 2. 

5. One-point correlation function and properties of the vacuum 

Let us consider the one-point correlation function 

(J(£,77)> = (0|exp(-^) J(fl expfaW) |0> = (0|J(0|0>. 

where J(£, ?7)-is conservation current in the theory. In a space-invariant 
theory, any one-point correlation function is equal to a constant 77)) = 
C. This property means that 

(0|J n |0) =5 nfi C, n = 0,±l,±2... . 

And we have the property that (0| and J n |0) are orthogonal to each other 
for all n^O. We can obtain this property if we assume that 

J n \0) = 0, fiGN. 
7 



In our case we have 

J*|0)=0, neN; keN. (5.1) 

We should have the properties of invariance of vacuum as follows: 

V\0) = 0, (5.2) 

where V-is the operator of translation (V — \ f+™ : (d^x(0) 2 '■ ^0- We 
know that 

[Vj k n ] = {k + n + l)J k n _ v (5.3) 

Now it is easy to obtain 

j fc |o) = o, JMo) = o, j* 2 |o) = o,... JVilo) = o, 

and for the left vacuum: 

(0|J fe = 0, (0|jf = 0, (0|J 2 fc = 0,... (0|J fc fe +1 = 0. 
It is easy to understand that 

(0\A 1 ^ l , Vl )....A n (UVn)J 1 p \0) =0, p = 0, ±1, ±2. 

We should know the commutator [J*, Ak(£k, Vk)]- We can calculate this 
commutator if A(£) = : exp(ax(0) : 

[J*, : exp(«x(0) ; ] = 

= C + V 1 -2,- expax(0 :)-{^ +1 (n+2)^+ ^ ~ a) (n+2)(n+l)H = expa*(0 : ■ 
Now we write the differential equations for the correlation function 

n 

^^ +2 (0|ii(^i,77i)...exp(-i?7 j 7Y) (Ji 2 ,: expa^xfe) exp^-W) ...A n (£„, r?„)|0) = 
j'=i 

= X>J/, p = 0,±l,±2, 

where are known differential operators. For the two-point correlation 
function, we have three equations 

( a 3 _ a _i_ 03 _ g)f 



«(6 - 6)(26 + 6)%/ + - 6X6 + 26)%/ 

-((a 3 -a)6 + (/3 3 -/3)6)/ = 0, 

«(6 - 6)(36 2 + 266 + 6 2 )%/ + mi - + 36 2 + 266)%/ 

-(2(a 3 -a)6 2 + 2(/3 3 -/3)6 2 )/ = 0. 
They have a solution (if we set (3 = —a) 

/(6,6^i,^) = C(a)G(r ?1 ,r ?2 )(6-6) (1 ^ 2) , 

where G(r]i,r] 2 ),C(a) are arbitrary functions. The equation (2.3) in light- 
cone coordinates has the symmetry (£ «-> 77), and so we have /(r) ~ y 2 ^-" ) 
where r is the distance between the points. We can investigate other Ward 
identities from another current J 2 (£) and obtain the same formula for f(r). 

Similar conservation currents exist in other models (see the Appendix). 
Now let us consider integrals of motion in systems with Hamiltonian 

/'+00 r+00 
■ exp(ax(0) : d£ + / : exp(fcax(0) : d f ■ 
-00 J — 00 

We have obtained nontrivial integrals of motion if k — — 1, — 2, , - 2 o, A and 
local conservation currents (see Appendix). When k = — 1 we have the sine- 
Gordon theory and if k — — 2, the Zhiber-Shabat theory. The theories with 
other k have no classical limit (a — > 0). 

6. Conclusion 

Computer calculations confirm an absolutely general behavior of quantum 
integrable systems with exponential interactions. We have obtained some 
examples of local conservated currents (densities) in all known models. Correlation 
functions in these cases can probably be obtained by a similar procedure. We 
have obtained a new infinite-dimensional symmetry in quantum integrable 
theories, and we calculated two-point correlation function using this symmetry. 

There is an interesting connection between this calculation and a similar 
calculation in fermionic models [6], where we used the Hamiltonian method to 
solve the massless Thirring model and obtained an anomaly for conservation 
of the energy-momentum tensor. We in fact reveal the hidden problem 
in fermionic anticommutators by demonstrating on obvious problem- the 
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anomaly. So the massless Thirring model is not a conformal field theory 
(solvable theory) and we can not use Coleman's way of proving the equivalence 
of the massive Thirring model and sine-Gordon theories using perturbative 
methods. 
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Appendix A. 

Here we report the arrangement and quantities m n of nontrivial currents 
at levels n. A current which is a differential of another current is trivial. The 
results of computer calculations are listed in the Table below. We investigate 
the space of quantum integrals of motion for the Hamiltonian Ti k . This 
problem can move to investigation of system S n ordinary linear equations 
with parameters a, k, degeneration of this system can correspond to appear 
the current at level n. Determinant of the system for n = 4, 6 have a form 

det S 4 = {a 2 - 1)(1 - k) 2 {a 2 k + 2) 2 {k + l){-a 2 k + 2), 

det Sq = 

(a - l) 2 (a + 1) 2 (3« 2 - 2)(a 2 - 3)(1 - k)\ka 2 + 2)\k + l)(Jfe + 2)x 

x (2k + l)(-ka 2 + 2)(-ka 2 + A)(-ka 2 + l)(-2a 2 k 2 + ka 2 - 2k - 3) 

After we find solutions for current when det S n = det S4 — O, which demand 
some conditions for parameters a, k. 

Table 
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a 


k 


m 2 


777.3 


m 4 


7775 


m 6 


m 7 


m 8 


7779 


77710 


777 n -conjecture 


1 


-1 





1 


1 


1 


1 


1 


1 


1 


1 


777„ = 1, for 77 > 2 


V2 


-1 


1 





1 





1 





1 





1 


"m 2 p — 1, p G N 


V3 


-2 














1 











1 


777 6+ 4 { p_i), pEM 




-2 














1 





1 





1 


777 6+ 2(p-l), P E N 




2 

« 2 




















1 





1 


to 8+2 (p-i), p e N 


Vs 


1 




















1 





1 


TO6+4(p-i), peN 


[2 

V 3 


1 

5* 














1 





1 





1 


777 6+ 2(p-l), p G N 
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